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SUMMARY 


The accurate prediction of loads on flexible, low aspect-ratio wings is critical to the 
design of reliable and efficient aircraft. The conditions for structural design frequently 
involve nonlinear aerod 3 mamics. 

Under previous NASA contracts (NASl-12875, NASl-14141, and NASI- 14962) a large 
experimental data base for three wing shapes was obtained, and linear theoretical 
methods were evaluated. The current contract, NASl-15678, extends the evaluation of 
state-of-the-art theoretical predictive methods to two separated-flow computer programs 
and also evaluates a semi-empirical method for incorporating the experimentally 
measured separated-flow effects into a linear aeroelastic analysis. 

The resultant three tasks have been documented separately. This volume describes the 
evaluation of R. P. White’s (RASA Division of Systems Research Laboratories) 
separated-flow method (Task I). The evaluation of The Boeing Company’s 
Three-Dimensional Leading-Edge Vortex (LEV) code (Task III) is presented in NASA 
CR-3642. The development and evaluation of a semi-empirical method to predict 
pressure distributions on a deformed wing by using an experimental data base (Task II) 
is described in NASA CR-3641. 

The method of R. P. White was developed for moderately swept wings with multiple, 
constant-strength vortex systems. The flow on the highly swept wing used in this 
evaluation is characterized by a single vortex system of continuously varying strength. 
The data comparisons, as currently formulated, show that this method does not predict 
the pressure distribution on this highly swept wing. 


INTRODUCTION 


Accurate anal 3 rtical techniques for the prediction of the magnitude and distribution of 
aeroelastic loads are required in order to achieve an optimum design of the structure of 
large flexible aircraft. Uncertainties in the characteristics of loads may result in an 
improper accounting for aeroelastic effects, leading to understrength or overweight 
designs and unacceptable fatigue life. In addition, the correct prediction of load 
distribution and the resultant structural deformation is essential to the determination 
of the aircraft stability and control characteristics, control power requirements, and 
flutter boundaries. The alternative to using satisfactory anal 3 rtical techniques is the 
increased use of expensive, time-consuming wind tunnel testing for each aircraft 
configuration. 

The problem of accurate load prediction becomes particularly acute for aircraft with low 
aspect-ratio wings where critical design conditions occur in the transonic speed regime. 
In this region, at t 3 q)ical design angles of attack, the flow is generally nonlinear - 
mixed flow, embedded shocks, separation, and vortex flow. 

A program was started in 1974 to systematically obtain experimental pressure data for 
an arrow wing throughout the subsonic, transonic, and low supersonic Mach numbers. 
This program was comprised of three NASA contracts: NASl-12875, NASl-14141, and 
NASl-14962 (documented in refs. 1 through 12). As the specific objective was to 
understand the change in load with aeroelastic deformation, three wing shapes were 
tested - all with the same planform and thickness distribution. The first wing was flat 
(no camber or twist); the second has a spanwise twist (typical of aeroelastic 
deformation) but no camber; and the third has the same twist with camber 
superimposed. 

In addition to the creation of a data base, which is useful for evaluating aeroelastic 
effects, a second objective was to evaluate state-of-the-art theoretical methods that 
might be used for this purpose. Primarily these methods were linear. The evaluations 
showed that linear theories are adequate at low angles of attack typical of cruise 
conditions and are basically capable of predicting loading changes due to smooth 
changes in wing shape at these low angles. However, at the higher angles of attack 
typical of structural design conditions, these methods are not useful because the flow is 
nonlinear due to leading-edge separation of the flow. The limited comparisons that were 
made with advanced separated-flow methods indicated some hope even though the 
aerodynamic panel model available at that time was very crude (only a few panels to 
represent the camber surface). 

The current evaluation of methods for predicting pressure distributions when the flow is 
separated is divided into three tasks. Two currently available computer codes were 
evaluated in Tasks I and III, and an approach involving semi-empirical corrections to 
linear theory was investigated in Task II. The three tasks are essentially independent 
efforts and are documented separately: Task I, an evaluaton of R. P. White’s computer 
code in this document; Task II, the development and evaluation of a semi-empirical 
method in NASA CR-3641; and Task III, an evaluation of Boeing’s Three-Dimensional 
Leading-Edge Vortex computer code in NASA CR-3642. 
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The computer code evaluated in Task I is the method of R. P. White of the RASA 
Division of Systems Research Laboratories, which models the vortex as a concentrated 
region of vorticity over the wing, and adds the results of this phenomenon to that of the 
linearized flow over a thin wing. The NASA arrow-wing configuration data base 
described above was used for this evaluation to determine the possible applicability of 
this method to the prediction of aeroelastic loads on highly swept wings. The following 
are discussed: theoretical basis of White’s method; Implementation of the method in a 
computer code; modeling of the arrow wing for analysis by the method; comparisons of 
calculated and experimental data; and the results of the evaluation. 
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SYMBOLS 


b 

BL 

c 

c, M.A.C. 

Cb 

Cc 

Cc 

Cm 

Cm 

Cm. 25c 

Cn 

Cn 

Cp 

D 

M 

MS 


Ps 


Pt 

q 

s 

Sh 


wing span, cm 

buttock line, cm; distance outboard from model plane of symmetry 

section chord length, cm 

mean aerodynamic chord length, cm 

surface bending moment coefficient referenced to y^^f, positive wingtip up 
surface chord force coefficient; positive aft 
section chord force coefficient; positive aft 

surface pitching moment coefficient, referenced to 0.25 M.A.C. ; positive 
leading edge up 

section pitching moment coefficient referenced to section leading edge; 
positive leading edge up 

section pitching moment coefficient referenced to section 0.25c; positive 
leading edge up 

surface normal force coefficient; positive up 

section normal force coefficient; positive up 

™ ^ measured pressure - reference pressure 

pressure coefficient = ^ 

body diameter, cm 
Mach number 

model station, cm; measured aft along the body centerline from the nose 
static pressure, kN/m^ 

total pressure, kN/m^ 
dynamic pressure, kN/m^ 

reference area used for surface coefficients, cm^ 

area of streamwise strip associated with a pressure station, cm^; used in 
summation of section force coefficients (app. B) 

free stream velocity 
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x,y,z general coordinates for distances in the longitudinal, lateral, and vertical 

directions respectively 

yj.gf distance outboard of model centerline of the bending moment reference 

point, cm 

a corrected angle of attack, degrees; the angle between the wing-root chord 

and the relative wind measured in the model plane of symmetry; includes 
compensation for sting deflection, tunnel-flow angularities, and wall 
effects; positive nose up with respect to relative wind 

“sec wing twist angle relative to wing reference plane, degrees; positive 

leading edge up 

/3 angle of sideslip, degrees; positive nose left with respect to relative wind 

S control surface deflection, degrees; positive leading edge down for leading 

edge (see exception in app. B) and trailing edge down for trailing edge 

Tj fraction of wing semispan, y/(b/2) 

A sweep angle, degrees; measured from a line perpendicular to the model 

centerline, positive aft 

(f> angle defining location of pressure orifices on the surface of the 

cylindrical body at a constant MS, degrees; measured from the top of the 
body 


Subscripts: 

L.E. leading-edge control surface 

r wing root 

s referenced to segment of local chord 

T.E. trailing-edge control surface 
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DATA BASE 


The data obtained, both experimental and theoretical, have been presented in several 
papers (refs. 1 through 3) and are presented in more detail in numerous NASA reports 
(refs. 4 through 12). 


WIND TUNNEL MODELS 

The configuration chosen for this study was a thin, low aspect-ratio, highly swept wing 
mounted below the centerline of a high fineness-ratio body. The general arrangement 
and characteristics of the model are shown in figure 1. Two complete wings were 
constructed for contract NASl-12875, one with no camber or twist and one with no 
camber but with a spanwise twist variation. A third wing with camber and twist was 
constructed for cbntract NASl-14962. Deflectable control surfaces were available on all 
three of these wings. 

The three wings, body, and fin used to create this data base are described in detail in 
appendix A. The wings all have the same planform, thickness distribution, and 
placement of orifices. The twisted wing and the cambered-tvristed wing have the same 
twist, i.e., the coordinates of the leading edges and trailing edges of the two wings are 
the same. This twist distribution is shown in figure 2. Sections at the root, midspan, and 
tip of the cambered-twisted wing (fig. 3) show not only the camber but the position of 
the sections of the cambered-twisted wing and the twisted wing, relative to the wing 
reference plane (flat wing). The flat wing had a sharp leading-edge segment in addition 
to the rounded leading-edge segment common to all three wings. 

The capability to measure the detailed load distribution on the wing and body of this 
configuration was provided by distributing 300 pressure orifices on the model. Each 
wing had 217 pressure orifices equally divided into seven streamwise sections on the 
left half. Orifices were located on both the top and bottom surfaces at the chord wise 
locations shown in figure 4. Pressure orifices were located on the body in five 
streamwise rows of 15 orifices each. An additional eight orifices in the area of the 
wing-body junction made a total of 83 orifices on the left side of the body. 

WIND TUNNEL TESTING 

The experimental data used in this study were obtained in the Boeing Transonic Wind 
Tunnel (BTWT) under NASA contracts NASl-12875 and NAS I- 14962. A description of 
the tunnel and tests are in appendix A. The current study was limited to the wings that 
had both leading-edge and trailing-edge control surfaces undeflected. Table 1 shows a 
summary of these data. 


DATA 

The measured pressures were edited, as necessary, to account for plugged or leaking 
orifices or missing data points. The pressure coefficients were then integrated, as 
described in appendix B, to obtain streamwise section coefficients and total surface 
coefficients. When pressure coefficients were required at points other than where 
measured, a linear interpolation was used. 
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Table 1. -Summary of Subsonic /Transonic Test Conditions by Test and Run Number 


Contract 

Test number 

Mach number 

Data 

document 

0.40 

I 

0.70 

0.85 

0.95 

1.00 

1.05 

1.10 

Flat wing, rounded leading edge 

NAS1 -12875 

BTWT 1415 

269 

263 

267 

266 

268 

264 

262 

NASA CR- 132727 

Flat wing, sharp leading edge 

NASI -12875 

BTWT 1415 

368 

366 

372 

374 

373 

367 

365 

NASA CR- 132727 


Twisted wing, rounded leading edge 


NASI -12875 
NASI -14962 

BTWT 1415 
BTWT 1627 

450 

15 

445 

14 

449 

13 

447 

12 

448 

11 

446 

10 

444 

NASA CR-1 32727 
NASA CR-1 65701 

Cambered-twisted wing, rounded leading edge (fin off) 

NASI -14962 

BTWT 1627 

43 

41 

40 

39 

38 

37 

- 

NASA CR-1 65701 




















































THEORY AND ASSUMPTIONS 


The basis of R. P. White’s prediction method (ref. 13) is that the flow over a wing at 
high angles of attack may be considered as the sum of the contributions of two quite 
separate phenomena: 

1) the flow over a thin wing, which is described by the linearized small disturbance 
potential equation, and 

2) the flow due to a concentrated region of vorticity, which occurs above the upper 
wing surface and originates due to the flow separation along the leading edge. 

These two phenomena are inconsistent as the presence of a region of vorticity violates 
the irrotational assumption of potential flow. 

Earlier investigators, however, have found that this assumption can be used with 
reasonable engineering accuracy. Brown and Michael (ref. 14), for example, represented 
the flow over a slender delta wing as a two-dimensional potential flow in the cross-flow 
plane, with a concentrated potential vortex added into the field above the wing. In their 
method, the strength and position of the vortex was determined by solution of the 
potential-flow equations subject to the boundary conditions of no flow through the 
surface, smooth flow at the leading edge (Kutta condition), and a condition of zero force 
on both the vortex and the feeding sheet transmitting vorticity from the wing leading 
edge to the concentrated vortex. Later investigators used similar methods with minor 
variations. A highly sophisticated version of this approach is the Boeing-developed 
Leading-Edge Vortex (LEV) method of reference 15, which uses a three-dimensional 
panel method for the wing surface and the rolled-up leading-edge vortex (evaluation 
discussed in Task III documentation, NASA CR-3642). 

The White method differs from these previous approaches in several aspects. Some of 
these may be described as follows: 

1) The position of the vortex is determined by the user and given as input data. This 
information may be determined from theoretical or empirical methods, such as that 
of Smith as described in reference 16, or from experimental data. The results of the 
program are very sensitive to the vortex location chosen. 

2) The strength of the vortex is determined iteratively from calculated leading-edge 
vortex strength given by an incompressible doublet-lattice panel method, not as 
part of a simultaneous solution of the coupled lifting-surface potential vortex flow 
as in the Brown and Michael method (ref. 14). The strength of the vortex is 
assumed constant in White’s method. 

3) The region of vorticity is considered to be of finite extent and is given by a 
self-similar solution of the Navier-Stokes equation. The suction pressures on the 
surface determined by the vortex solution are superimposed on those calculated by 
the doublet-lattice lifting-surface method. 

4) The vortex is considered to be free, i.e., not connected to the wing surface by a 
feeding sheet transmitting vorticity from the leading edge to the vortex core. Thus, 
the Kutta condition of smooth flow off the leading edge cannot be inherently 
satisfied. 



Further details of White’s method, including a description of the viscous flow solutions 
and the empirical constants used, may be found in reference 13. 

The method described in reference 13 has been incorporated in a computer code for use 
on the Control Data Corporation 6600/Cyber 175 series of computers. This code is 
described fully in reference 17. During the performance of the contract, the code was 
received from RASA, installed on the Boeing computers, checked out and validated, and 
modified as required for the NASA arrow-wing configuration. 

VALIDATION OF THE COMPUTER CODE 

The case used by RASA for development and validation of the code was a moderate 
aspect-ratio, moderately swept planform with a strake based on the F-4 airplane (see 
ref. 13). The flow over this planform at high angles of attack is characterized by the 
presence of multiple vortex systems (strake, leading edge, and wing tip) that tend to 
increase in strength only over a small portion of the lifting surface and behave as free 
vortices over most of the wing. An excellent comparison of measured and calculated 
pressures was obtained by RASA for the F-4 wing. 

The F-4 planform was used as a test case to validate the code after its installation at 
Boeing. After receiving additional corrections from RASA, the calculated results were in 
agreement with those provided by RASA (within the accuracy of computer and 
operating system differences). 

MODIFICATION OF THE COMPUTER CODE 

The computer code, as received, needed several modifications before application to the 
NASA arrow wing could be attempted. These modifications were made jointly by Boeing 
and RASA personnel and may be described as follows: 

1) The airfoil section definition in the code, as received from RASA, is given by a 
polynomial definition (see ref. 18) representing the NACA OOXX series of airfoils 
(where XX is the thickness ratio in percent chord). This airfoil is assumed constant 
across the span. Different airfoil definitions were incorporated in the code by 
including the polynomial coefficients obtained from a least squares fit of the 
surface ordinates used to construct the arrow-wing wind tunnel models. 

2) In the original code, the mean plane of the wing was assumed to be flat. For the 
twisted and cambered-twisted arrow-wing configurations, modification of the 
surface geometry to account for spanwise variations was required. 

3) As received, the code did not contain an iteration procedure to automatically 
determine vortex strength from the potential-flow calculation. This was 
incorporated by RASA during the contract. 

Other modifications, which were discussed but not made, were the incorporation of 
compressibility (the current code is for incompressible flow), inclusion of a feeding sheet 
model, iteration of the vortex position as well as strength (this was included in the code 
but not found to be operative), and modifications to the separated flow criteria described 
in reference 13. 
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APPLICATION OF THE METHOD 


After installation and validation of the code, the NASA arrow- wing configuration was 
modeled for input to the program. The cases that were run for comparison with the 
experiment consisted of the flat, twisted, and cambered-twisted wings at a Mach 
number of 0.40 and an angle of attack of 16®. In addition, the flat wing was run at an 
angle of attack of 8®. The results of these cases and the comparison with experimental 
data are described in this section. 


MODELING 

The modeling consists of three general parts. The procedure for this modeling, along 
with the difficulties that were encountered, are described. 

PLANFORM GEOMETRY AND PANELING 

Surface paneling of the wing is required for implementation of the doublet-lattice 
method for the linear part of the flow. Paneling in the code is either manual (user input 
of constant percent chord (x/c) and constant fraction of semispan (y/(b/2)) locations for 
the panel edges) or automatic (a standard set of nondimensional paneling). The 
automatic option was used for the arrow wing. Figure 5 shows the paneling used. The 
code assumes the paneling is the same for the upper and lower surfaces of the wing. 

SURFACE GEOMETRY 

In the White method, the boundary condition of no flow through the surface is specified 
on the actual airfoil surface. It is, therefore, necessary to provide an accurate 
representation of the surface thickness, camber, and twist distributions. In the original 
code, an NACA OOXX series airfoil was used with the surface shape represented by the 
polynomial description given in reference 18. The wing was assumed to be symmetric 
about the mean plane, i.e., no twist or camber. In order to represent the various 
configurations of the NASA arrow wing, development of new surface representations 
and modifications to the code were required. 

The polynominal coefficients to represent the thickness of these airfoils (with both 
rounded and sharp leading edges) were obtained with a least squares fit of the surface 
coordinates of the airfoil (z/c versus x/c). Table 2 shows the form of the polynomial used 
to fit the geometry and the coefficients obtained. These coefficients were used to define 
surface geometry and slope in White’s code. A comparison of the forward portions of the 
rounded- and sharp-airfoil sections, as used, are shown in figure 4. 

The nondimensional airfoil sections also varied as a function of span because the 
trailing edge is a constant 0.0254 cm (0.01 in.) thick due to manufacturing constraints 
on the models. It was assumed for purposes of incorporating the geometry in White’s 
code that a constant spanwise section could be used. This assumption caused some 
difficulty as the answers varied greatly, depending on whether the root section or tip 
section was chosen as representative. 

The root section, which has a relatively sharp trailing edge, was originally used to 
define the airfoil. This gave large bumps in the pressure distribution near the trailing 
edge, which were not observed in the experimental data. As previously mentioned, the 
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Table 2.-Polynomial Coefficients for Airfoil Definition 



Ai + A2 


VF^ 


10 

^ An +2 
n=l 


(■!■)" (i) 




Rounded LE. 
root section 

Rounded L.E. 
tip section 

Sharp L.E. 
tip section 


Ai 

0.00 

0.00 

0.00 


A2 

.10095021 E +00 

.10093185E+00 

0.00 


A3 

-.20101451 E+00 

-.19951709E+00 

.26196052E+00 


A4 

.47539204E+00 

.45661 141 E+00 

-.10548058E+01 

Upper 

A 5 

-.22656037E+01 

-.20088269E+01 

-.54993659E+01 

surface 

Ab 

.1 144831 1E+02 

.97680669E+01 

.68051393E+02 


A7 

-.38620912E+02 

-.32431 182E+ 02 

-.28466430E+03 


Ae 

.81190343E+02 

.67393001 E+ 02 

.65260528E+03 


A& 

-.10711414E+03 

-.88096659E+02 

-.89728052E+03 


A10 

.86521 11 6E +02 

.70665296E+02 

.73847873E+03 


A11 

-.39099350E+02 

-.31768977E+02 

-.33591 049E+03 


Ai2 

.75650206E+01 

.61223848E+01 

.65013256E+02 


tip section had a slight bluntness at the trailing edge when measured in percent of 
chord because of manufacturing constraints on the model. Figure 6 shows the 
nondimensional geometry of the trailing edge of the root and tip sections with the 
vertical scale exaggerated. When the tip section (blunt trailing edge) was used, the 
bumps in the pressure distribution were not observed. A comparison of the resulting 
pressures for the two airfoil definitions is shown in figure 7. The cause of the 
discrepancy is not known at this time, but it is suspected to be anomalies in the 
numerical method for calculating the linear part of the pressure distribution. The 
modeling for the subsequent comparisons used the definition that was based on the tip 
airfoil. 

Camber and twist distributions were added to the basic symmetrical surface to obtain 
models for comparison with experimental data on the twisted and camber ed-twisted 
wings. These distributions for the model are defined in appendix A. 

VORTEX MODELING 

The vortex modeling requires both the vertical and horizontal location of the vortex. In 
the analysis of the F-4 wing performed by RASA, flow visualization techniques were 
used to determine the horizontal and vertical locations. Three-dimensional 
flow-visualization data were not available for the arrow wing, although surface pressure 
data could be used to determine the horizontal location. Only a single vortex system 
was modeled for the arrow wing, whereas three systems were used for the F-4. 
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Several different approaches were tried in modeling the vertical position of the vortex. 
The initial approach was to obtain the vertical position by assuming that the angle 
between the wing and the vortex was 
one half the angle of attack (see 
sketch), i.e., 8® for an angle of attack of 
16®. This is the theoretical limit for a 
concentrated force-free vortex above a 
slender delta wing of vanishing aspect 
ratio. For this assumption, the vortex 
had almost no effect on the wing 
surface pressures, and considerable 
disagreement was observed in the 
comparison of calculated and measured 
pressures. 




The second approach was to use a vortex similar to that used for the F-4 wing; 
one half the angle of attack to the 0.25 root-chord position, then at a constant 
above the surface to the trailing edge, then parallel to the free stream aft of the 
edge (see sketch). With this modeling, 
a significant effect of the vortex on the 
surface pressures was observed. 

However, this effect appeared to be 
exaggerated compared with 
experimental pressures, i.e., more 
suction pressure due to the vortex. 


namely, 

distance 

trailing 




At the suggestion of RASA personnel, the third approach was to use the method 
described by Smith (ref. 16) to locate the vortex position. This method assumes conical 
flow and, therefore, provides a straight vortex emanating from the apex of the wing and 
includes the effect of both angle of attack and semi-apex angle. Starting the vortex at 
the wing apex gives a very large, unrealistic suction peak. Therefore, the vortex was 
started at a location 20 percent of the root chord aft of the apex. As shown in figure 8, 
at 16® angle of attack, the horizontal location of the vortex, as given by the Smith 
method, agreed excellently with that determined from the suction peaks in the 
experimental data for the flat wing with rounded leading edge. The vertical location of 
the vortex (Smith method) shows the vortex above the wing at approximately one 
quarter of the angle of attack. For 8® angle of attack, however, the pressure data 
indicates that the vortex starts along the leading edge at about the midspan point, and 
proceeds aft in a curving path to the trailing edge, with a fairly significant 
disagreement with the conical flow assumption, as shown in figure 8. 


The Smith method for placing the vortex was used for all calculated results. 


In White’s method, the leading edge vortex is assumed to contain a central core region 
of predominantly axial turbulent flow. In the core, the circumferential velocity is 
assumed to vary directly with vortex radius, resembling a solid-body rotation. Outside 
the core, the circumferential velocity varies inversely with radius, similar to a potential 
vortex. The core radius is given by the solution to the viscous flow equations, but an 
initial value must be chosen. The initial core radius is a user-input parameter and is 
ideally obtained from photographic flow-visualization data which, however, were not 
available for the arrow wing. A nominal value of 0.0107 m (0.035 ft) was initially 
chosen. In an attempt to improve the comparisons, a sensitivity study on the initial core 
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radius was performed. Higher and lower values of 0.0427 m (0.140 ft) and 0.0030 m 
(0.01 ft) were tried. A comparison of the calculated pressures for these three values is 
shown in figure 9. It may be seen that no consistent trend in the pressures can be 
observed due to this change in core radius, therefore, the initial core radius of 0.0107 m 
(0.035 ft) was used for the study. 


COMPARISON OF ANALYSIS AND TEST 

Calculations on the arrow-wing configuration were performed using White’s code, 
modified to account for the exact surface geometry. Since the code was limited to 
incompressible flow, the results are compared to experimental data at a Mach number 
of 0.40, the lowest Mach number tested. Most of the calculations were done for an angle 
of attack of 16®. Some calculations were later performed at about 8®. 

Chordwise pressure distributions on the arrow wing were measured at seven spanwise 
locations: y/(b/2) = 0.09, 0.20, 0.35, 0.50, 0.65, 0.80, and 0.93. The code, as written, 
evaluates pressures at y/(b/2) = 0.01 and at each increment of 0.02 across the wing. This 
difference in fraction of semispan of 0.01, where it occurs, was assumed negligible in 
comparing calculated and experimental pressures. 

Figures 10 through 13 show a comparison of calculated and experimental data at 16® 
angle of attack. Figure 10 shows the calculated and measured upper-surface pressure 
coefficients for the flat wing with sharp leading edge. The pressure distribution is not 
predicted accurately and is not adequate for design loads. The calculated pressure 
distributions show "humps,” which may be interpreted as the effect of vortex suction. At 
one section (y/(b/2) = 0.20), the comparison of theoretical and experimental pressures 
shows fairly good agreement; unfortunately, this appears to be merely fortuitous. The 
vortex effect is far more pronounced in the experimental data, particularly at midspan 
(y/(b/2) = 0.35 and 0.50), where the values of the theoretical data are much lower than 
the experimental. On the outboard wing (y/(b/2) = 0.65, 0.80, and 0.93), the 
experimental data indicate that the flow is separated, i.e., the pressure level is nearly 
constant without significant pressure recovery near the trailing edge. The calculated 
values, on the other hand, are typical of classical linear theory, with high suction 
pressures near the leading edge and recovery to static pressure near the trailing edge. 

White’s method contains a "separated-flow criteria” for determining whether, at any 
point on the wing, separation has occurred. A section was assumed to be stalled if the 
net aerodynamic angle of attack (the geometric angle of attack minus the induced 
angle) was greater than an empirically defined angle. Considerable effort on this 
separated flow criteria failed to reproduce the experimentally observed behavior. It 
would appear, therefore, that significant difficulties in modeling the vortex flow remain 
in the code. Similar observations can be made for the flat, twisted, and 
cambered-twisted wings with the rounded leading edge, shown in figures 11 through 13. 

A comparison of experimental and calculated data for the flat wing with rounded 
leading edge at 8® angle of attack is shown in figure 14. The comparison at 8® is better 
than that at 16®, probably because at 8® the vortex is present only over the outboard 
portion of the wing (fig. 8), and it is weaker as well. This is consistent with the observed 
difficulties in predicting the effect of the vortex. Since at 8® angle of attack the flow is 
predominately potential, the vortex would be expected to have a smaller effect on the 
lifting-surface pressure relative to the linear portion predicted by the doublet-lattice 
method at 8® than at 16® angle of attack. 
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CONCLUSIONS 


White’s method has been evaluated for its applicability to the calculation of surface 
pressures that are due to separated vortex flow on a highly swept configuration using a 
variety of surface shapes. The shapes are those that would be t 3 ^ical of the deformed 
geometry of a low aspect-ratio wing under aeroelastic load. 

From the comparisons of experimental and calculated surface pressures, it does not 
appear that the method is capable of accurate predictions of pressure with the modeling 
that is available in the code. The method was developed for applicability to less highly 
swept wings with multiple vortex systems. These vortices were not characterized by the 
continuous variation of vortex strength that is typical of highly swept configurations. 
Furthermore, the method is highly sensitive to the placement of the vortex in the flow 
field. Although this information could be made available from wind tunnel tests for an 
undeformed wing, the use of the method in an aeroelastic solution would require that 
the solution be capable of predicting changes in the vortex position in response to 
surface geometric changes due to aeroelastic deformation. 

It cannot be recommended, therefore, that the method be considered further for 
incorporation in an aeroelastic analysis unless significant changes are made to: 

o the assumptions used in modeling the vortex system, especially to allow changes in 
vortex strength, 

0 the solution procedure to allow the vortex to be properly repositioned automatically 
as part of the solution, and 

0 include the effects of compressibility. 

Boeing Commercial Airplane Company 
P. O. Box 3707 

Seattle, Washington 98124 
May 1982 
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APPENDIX A 


DESCRIPTION OF DATA BASE 


WIND TUNNEL MODELS 

The configuration chosen for this study was a thin, low aspect-ratio, highly swept wing 
mounted below the centerline of a high fineness-ratio body. The general arrangement 
and characteristics of the model are shown in figure A:-l. Two complete wings were 
constructed for contract NASl-12875, one with no camber or twist and one with no 
camber but with a spanwise twist variation. A third wing with camber and twist was 
constructed for contract NASl-14962. Deflectable control surfaces were available on all 
three of these wings. 

FLAT WING 

The mean surface of the flat wing is the wing reference plane. The nondimensional wing 
thickness distributions (shown in table A-1) deviate slightly from a constant for all 
streamwise sections to satisfy a manufacturing requirement for a finite thickness of 
0.0254 cm (0.01 in.) at the trailing edge. The wing was designed with a full-span, 
25-percent chord, trailing-edge control surface. Sets of fixed angle brackets allowed 
streamwise deflections of ±4.1®, ±8.3®, ±17.7®, and ±30.2®, as well as 0.0®. A removable 
full-span leading-edge control surface (15 percent of streamwise chord) could be placed 
in an undeflected position and also drooped 5.1® and 12.8® with fixed angle brackets. 
Both the leading- and trailing-edge control surfaces extended from the side of the body 
(0.087 b/2) to the wingtip and were split near midspan (0.570 b/2). The inboard and 
outboard portions of the control surfaces were able to be deflected separately and were 
rotated about points in the wing reference plane. An additional leading-edge control 
surface for this wing was constructed with a sharp (20® included angle) leading edge to 
examine the effects of leading-edge shape. The surface ordinates and slopes of this 
leading-edge segment were continuous with those of the flat wing at the leading-edge 
hingeline (table A-1). The sharp leading edge was smoothly faired from 0.180 b/2 into 
the fixed portion of the rounded leading edge at 0.090 b/2. 

TWISTED WING 

The mean surface of the twisted wing was generated by rotating the streamwise section 
chord lines about the 75-percent local chord points (trailing-edge control surface 
hingeline). The spanwise variation of twist is shown in figure A-2. The hingeline was 
straight and located in the wing reference plane at its inboard end (0.087 b/2) and 
2.261 cm (0.890 in.) above the wing reference plane at the wingtip. The airfoil thickness 
distribution (table A-1) and the trailing-edge control surface location and available 
deflections were identical to those of the flat wing. 

CAMBERED-TWISTED WING 

The mean surface of the cambered-twisted wing was generated by superimposing a 
camber on the twisted-wing definition but keeping the coordinates of the leading edge 
and trailing edge of the cambered-twisted wing the same as those of the twisted wing. 
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Table A- l.—Wing Half -Thickness Distribution, Percent Chord 
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The camber is defined analytically in two parts: a typical cruise airfoil (basic) camber; 
and, an estimate of the aeroelastic deformation at a moderate positive angle of attack. 
The aeroelastic deformation was based on calculations - using a typical configuration - 
of deformation under load. This definition was modified slightly to provide zero camber 
at the model centerline so this wing would fit on the existing model body. This was 
achieved by using a factor k (fig. A-3) on the basic camber term, which provides a 
transition from no camber at the model centerline to the definition camber at 0.25 b/2. 
The defined wing is smoother than indicated in this figure as section geometry was 
directly calculated at only those sections marked in figure A-3. The full equation for the 
camber is: 



The resulting nondimensional camber is shown in table A- 2. The camber at the tip is 
approximately a 6® arc of a circle with the leading and trailing edges up. Sections at the 
root, midspan, and tip (fig. A-3) show not only the camber but the position of the 
sections of the cambered- twisted wing and the twisted wing relative to the wing 
reference plane (flat wing). The airfoil thickness distribution (table A-1) and the 
trailing-edge control surface location and available deflections were identical to those of 
the flat wing. 

BODY 

The body was circular in cross section and had a straight centerline. The body geometry 
is shown in figure A-1. The sting was an integral part of the model body. 


RELATIVE WING AND BODY LOCATION 

The wing reference plane was located 3.149 cm (1.240 in.) below and parallel to the body 
centerline (zero incidence). The apex (extension of the wing leading edge to the 
centerline) of the wing was located 33.496 cm (13.187 in.) aft of the model nose. 
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WING FIN 


The wing fin is a 3 -percent biconvex airfoil placed streamwise and perpendicular to the 
wing reference plane on the upper surface of the cambered-twisted wing at 0.725 
semispan. The dimensions of the fin and its relationship to the wing are shown in 
figure A-4. To obtain configurations with the outboard tr ailing-edge control surface 
deflected, the fin is extended down to touch the top of the control surface. (See fig. A-4.) 

PRESSURE ORIFICE LOCATIONS 

All pressure orifices were located on the left side of the model and distributed as shown 
in figure A-5 and tables A-3 and A-4. The flat wing with rounded leading edge, the 
twisted wing, and the cambered-twisted wing each had 214 orifices distributed in 
streamwise pressure stations of 31 (or 30) orifices at each of seven spanwise locations. 
One of these orifices was located at the leading edge; the remainder were distributed so 
that upper- and lower-surface orifices were located at the same chordwise locations. The 
orifice locations on the sharp leading edge were identical except for the omission of the 
leading-edge orifice at each spanwise station. The 83 orifices on the body were located 
at 15 stations along the length of the model. At each station, orifices were located at 
angles of 0®, 45®, 90®, 135®, and 180® measured from the top of the body. In the area of 
the wing-body intersection, the orifices that are nominally identified as being at 135® 
and 180® were located on the wing lower surface at the same lateral location as the 
orifices at 45® and 0®, respectively. Eight additional orifices were placed on the body, 
close to the juncture of the body with the wing upper surface. 

DESIGN AND CONSTRUCTION 

The objectives of this study dictated that the contours and physical characteristics of 
the flat wing, the twisted wing, and the cambered-twisted wing be as nearly identical as 
possible. The model was constructed of steel to minimize aeroelastic deflections and to 
provide strength for testing to a Mach number of 3.0. The aft body was flared 
approximately 4° from 194.310 cm (76.500 in.) aft of the nose to provide the required 
safety factor on predicted loads (fig. A-1). The model size was selected as the best 
compromise between minimizing potential tunnel blockage and providing adequate 
room to install orifices in the model. 

A computerized lofting program was used to provide the wing definition. This definition 
was then used to machine the model components using numerically controlled 
machines. The tolerance on the contour was +0.1524, -0.0 mm (+0.006, —0.0 in.). The 
leading- and trailing-edge control surfaces were cut from the wings after they had been 
machined to final contour. A cut along the 15-percent chord line of the twisted wing 
removed enough material to simulate the elastic characteristics of the flat wing 
(fig. A-6). As a result of the previous tests it was determined that it was not necessary 
to remove this material on the cambered-twisted wing as the wings were very rigid. 
Fixed angle brackets (arranged as shown in fig. A-6) were used to obtain the required 
control surface deflections with all pivot points located midway between the upper and 
lower surfaces at the hingelines. The brackets were also machined on numerically 
controlled machines. The same sets of trailing-edge brackets were used on all three 
wings, and the same sets of leading-edge brackets were used for both the rounded and 
sharp leading edges. 
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Table A-3.— Wing Pressure Orifice Locations, Percent Local Chord 


(a) Section at 0.09-^, chord = 102.89 cm 



Nominal 

Flat wing, 

a = 0.0° 
sec 

Twisted wing, 
«sec = -0'01° 

Cambered-twisted wing, 
“sec “-0-01° 

Rounded leading edge 

Sharp leading edge 

Rounded leading edge 

Rounded leading edge 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

0.00 

0.00 




— 

0.00 


0,00 


2.50 

2.45 

2.59 

2.61 

2.54 

2.26 

2.26 

2.58 

2.51 

5.00 

4.95 

5.07 

5.06 

5.03 

4.76 

4.76 

5.10 

5.04 

8.50 

8.45 

8.53 

8.59 

8.58 

8.40 

8.26 

8.64 

8.56 

11.30 

— 

— 

— 

11.31 

— 

— 

- 

— 

12.25 

— 

— 

— 

— 

12.23 

12.27 

— 

— 

12.50 

12.45 

12.55 

12.58 

— 

— 

— 

12.63 

12.54 

17.50 

17.49 

17.62 



17.59 

17.66 

17.64 

17.55 

20.00 

19.94 

20.08 



20.03 

20.03 

20.14 

20.00 

30.00 

29.92 




29.98 

29.89 

30.14 

30.00 

45.00 

45.00 




44.96 

44.89 

45.12 

45.03 

00.00 

59.98 




60.01 

59.97 

60.11 

60.00 

70.00 

70.03 

70.13 



70.05 

69.95 

70.09 

70.04 

72.50 

72.55 

72.60 



72.58 

72.51 

72.62 

72.54 

77.50 

77.53 

77.62 



77.56 

77.51 

77.63 

77.52 


85.11 

85.14 



85.03 

85.00 

85.12 

85.04 


90.10 

90.10 



90.04 

89.98 

90.12 

90.00 

1 95.00 

95.09 

95.05 



94.96 

94.98 

95.10 

95.03 
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Table A-3.—( Con tinued) 


(b) 


Section at 0.20~, chord = 91 .80 cm 
2 



Nominal 

Flat wing, 
^sec = 0-0° 

Twisted wing, 

«sec = -0.47° 

Cambered-twisted wing, 
^sec = -0-47° 

Rounded leading edge 

Sharp leading edge 

Rounded leading edge 

Rounded leading edge 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

■1 

0.00 




0.00 


0.00 



2.59 

2.69 

2.62 

2.65 

2.52 

2.42 

2.63 

2.59 


5.05 

5.00 

5.14 

5.14 

5.00 

4.93 

5.09 

5.05 

8.50 

8.54 

8.59 

8.67 

8.62 

8.52 

8.40 

8.61 

8.64 

11.40 

— 

— 

— 

11.37 

— 

— 

— 

— 

12.50 

12.54 

12.49 

12.63 

— 

12.53 

12.42 

12.51 

12.62 

17.50 

17.63 

17.61 



17.65 

17.52 

17.59 


20.00 

20.08 

20.07 



20.00 

19.90 

19.95 


30.00 

30.04 

30.09 



30.02 

29.89 

30.05 


45.00 

45.08 

45.09 



45.03 

44.92 

45.04 


60.00 

60.02 

60.13 



60.03 

59.91 

60.02 

60.06 

70.00 

70.11 

70.13 



70.06 

69.96 

70.03 

70.01 

72.50 

72.63 

72.61 



72.55 

72.50 

72.59 

72.67 

77.50 

77.59 

77.65 



77.59 

77.52 

77.53 

77.57 

85.00 

85.07 

85.13 



85.02 

85.00 

85.09 

85.10 


90.14 

90.11 



90.07 

89.97 

90.04 

89.98 

1 95.00 

95.14 

95.10 



95.05 

95.08 

95.06 

94.98 
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Table A-3.—( Con tinued) 


(c) 


Section at 0.35— chord = 76.69 cm 
2 



Nominal 

Flat wing, 
Sec = 0-0° 

Twisted wing. 
Sec =-1-70" 

Cambered-twisted wing. 
Sec = -1-70° 

Rounded leading edge 

Sharp leading edge 

Rounded leading edge 


Rounded leading edge 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 


0.00 


— 



0.00 


0.00 

■M 


2.45 

2.59 

2.59 

2.58 

2.39 

2.33 

2.76 



4.93 

5.07 

5.11 

5.04 

5.12 

4.78 

5.05 


8.50 

8.60 

8.54 

8,65 

8.63 

8.49 

8.32 

8.68 

8.70 

10.50 

— 

— 

— 

10.46 

— 

— 

— 

— 

11.00 

— 

11.03 

— 


— 

— 

— 

— 

12.50 

12.37 

— 

12.57 

— 

12.50 

12.33 

12.59 

12.68 

17.50 

17.64 

17.63 



17.54 

17.53 

17.64 

17.62 

20.00 

20,00 

20.09 



19.94 

19.84 

20.03 

20.07 

30.00 

30.01 

30.10 



29.88 

29.87 

30.00 

29.93 

45.00 

44.99 

45.09 



44.96 

44.79 

45.00 

45.13 

60.00 

60.03 

60.08 



59.97 

59.89 

60.00 

60.10 

70.00 

70.07 

70.08 



70.03 

69.90 

70.04 

70.03 

72.50 

72.55 

72.58 



72.56 

72.44 

72.61 

72.52 

m 


77.61 



77.54 

77.51 

77.50 

77.60 



85.14 



85.08 

84.96 

85.09 

84.93 

Bin 


90.09 



89.89 

89.89 

89.98 

90.04 

1 95.00 

95.07 

95.09 




94.95 

94.86 

94.98 

95.10 
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Table A-3.— (Continued) 


(d) Section at 0.50—, chord = 61.57 cm 
2 


Wing reference 
plane 


Wing reference 
plane 


Nominal 

Flat wing, 

a = 0.0° 
sec 

Twisted wing, 

“sec =-2-85° 

Cambered-twisted wing, 

“sec =-2.85° 

Rounded leading edge 

Sharp leading edge 

Rounded leading edge 

Rounded leading edge 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 


0.00 




_ 

0.00 


0.00 



2.47 

2.53 

2.69 

2.60 

2.44 

2.38 

2.78 

2.62 


4.99 

4.95 

5.13 

5.06 

4.92 

4.80 

5.13 

5.15 

8.50 

8.48 

8.38 

8.66 

8.61 

8.46 

8.38 

8.64 

8.56 

10.10 

— 

— 

— 

10.14 

— 

— 

— 

— 

11.10 

— 

11.08 

— 

— 

— 

— 

— 

— 

12.50 

12.39 

— 

12.61 

— 

12.50 

12.31 

12.71 

12.55 

17.50 

17.64 

17.52 



17.54 

17.24 

17.71 

17.44 

20.00 

19.98 

19.97 



19.92 

19.83 

20.15 

19.89 

30.00 

30.07 

30.06 



29.91 

29.85 

30.04 ■ 

29.72 

45.00 

44.98 

45.06 



45.00 

44.85 

44.95 

44.97 


59.97 

60.00 



59.95 

59.92 

59.96 

59.94 

70.00 


70,10 



70.03 

69.88 

69.93 

69.86 

72.50 

wsm 

72.61 



72,56 

72.44 

72,53 

72.34 

77.50 

77,66 

77.65 



77.61 

77.43 

77.58 

77.43 

85.00 

85.19 

85.18 



84.85 

84.90 

84.96 

84.92 

90.00 

90.22 

90.12 



89.93 

89.93 

89.94 

89.91 

95.00 

95.05 

94,94 



94.88 

94.93 

94.98 

94.88 












































Table A-3.— (Continued) 


(e) 


Section at 0.65—, chord = 46.46 cm 
2 




Flat wing, 

a = 0.0° 
“sec 

Twisted wing. 

Cambered-twisted wing, 

«sec=-3-59° 

Nominal 

Rounded leading edge 

Sharp leading edge' 

Rounded leading edge 

Rounded leading edge 


Upper 

Lower 

Upper 

Lower 

Upper 

Lower 

Upper 

Lower 


surface 

surface 

surface 

surface 

surface 

surface 

surface 

surface 

0.00 

0.00 




0.00 


0.00 


2.50 

2.56 

2.66 

2.49 

2.38 

2.18 

2.49 

2.76 

2.79 

5.00 

5.06 

5.12 

4.94 

4,95 

4.76 

5.01 

5.35 

5.40 

8.50 

8.55 

8.55 

8.46 

8.40 

8.32 

8.45 

8.85 

8.74 

12.20 

- 

— 

12.12 

— 

12.21 

— 

— 

— 

12.60 

12.57 

— 

— 

— 

— 

— 

12.71 

12.77 


17.60 

17.65 



17.24 

17.44 

17.74 

17.58 


20.17 

20.11 



19.70 

19.88 

20.19 

19.96 


30.05 

30.11 



30.26 

29.73 

30.13 

29.85 

45.00 

45.16 

45.23 



44.75 

44.89 

45.03 

44.75 

60.00 

60.13 

60.13 



59.81 

59.87 

60.02 

59.99 

70.00 

69.89 

70.12 



69.92 

69.90 

70.09 

69.88 

72.50 

72.59 

72,69 



72.38 

72.49 

72.83 

72.15 

77.50 

77.74 

77.76 



77.22 

77.49 

77.56 

77.43 

85.00 

85.25 

85.32 



84.79 

84.93 

84.93 

84.76 

90.00 

90.22 

90.21 



89.70 

89.92 

89.95 

89.98 

95.00 

95.13 

95.27 



95.12 

94.86 

94.97 

94.98 
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Table A-3.—(Con tinued) 


(f) 


Sectfon at 0.80—, chord = 31.35 cm 
2 



Nominal 

Flat wing, 

a = 0.0° 
sec 

Twisted wing, 
“sec= -3.84° 

Cambered-twisted wing, 

“sec = -3-84° 

Rounded leading edge 

Sharp leading edge 

Rounded leading edge 

Rounded leading edge 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

Upper 

surface 

Lower 

surface 

0.00 

0.00 




0.00 


0.00 


2.50 

2.55 

2.47 

2.50 

2.46 

2.33 

2.43 

2.76 

2.62 

5.00 

5.01 

5.02 

5.01 

4.93 

4.86 

4.74 

5.27 

5.21 

8.50 

8.55 

8.59 

8.58 

8.41 

8.32 

— 

8.78 

8.54 

12.50 

12.50 

— 

12.58 

— 

12.47 

12.43 

12.69 

12.58 

17.50 

17.53 

17.57 



17.36 

17.47 

17.83 

17.34 

20.00 

20.16 

20.13 



19.79 

19.82 

20.11 

19.79 

30.00 

30.00 

30.11 



29.83 

29.83 

30.15 

29.48 

45.00 

44.91 

45.15 



44.81 

44.91 

44.81 

44.75 

60.00 

59.94 

60.10 



59.80 

59.92 

59.84 

59.79 

70.00 

70.06 

70.11 



69.89 

69.87 

69.77 

69.94 

72.50 

72.61 

72.60 



72.22 

72.39 

72.50 

72.33 

77.50 

77.73 

77.72 



77.29 

77.41 

77.22 

77.40 

85.00 

85.25 

85.18 



84.80 

84.95 

84.92 

84.92 

90.00 

90.20 

90.34 



90.62 

90.03 

90.19 

90.09 

95.00 

95.41 

95.49 



95.71 

95.00 

95.05 

94.94 
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Table A-3.—( Concluded) 


(g) 


Section at 0.93-^, chord = 18.25 cm 


Wing reference 
plane 


T.E, hingeline 




Flat wing, 
®sec ” 

Twisted wing, 

a = _4 1 4° 
“sec 

Cambered-twisted wing, 
“sec = -4-14“ 

Nominal 

Rounded leading edge 

Sharp leading edge 

Rounded leading edge 

Rounded leading edge 


Upper 

Lower 

Upper 

Lower 

Upper 

Lower 

Upper 

Lower 


surface 

surface 

surface 

surface 

surface 

surface 

surface 

surface 

0.00 

0.00 








2.51 

1.70 

1.81 

2.12 

1.86 

1.74 

2.59 


2.26 

5.00 

4.38 

4.68 

4.72 

4.52 

4.41 

4.65 


4.79 

8.50 

7.89 

8.24 

8.21 

8.06 

7.92 

8.23 

8.64 

8.13 

11.59 

— 

— 

— 

— 

11.59 

— 

— 

— 

12.25 

12.33 

— 

12.19 

— 

— 

— 

12.64 

12.16 

17.50 

17.36 

16.60 



16.60 

17.49 

18.03 

16.83 

20.00 

19.78 

19.81 



19.58 

19.96 

19.94 

19.44 

30.00 

29.67 

29.00 



29.17 

29.62 

30.22 

28.66 

45.00 

44.70 

44.80 



44.12 

44.44 

44.33 

44.77 

60.00 

59.68 

59.47 



59.18 

59.71 

59.47 

59.38 

70.00 

69.69 

70.33 



68.99 

69.31 

69.10 

70.07 

72.50 

72.15 

71.89 



71.59 

72.01 

71.78 

72.74 

nsi 

77.38 

77.31 



76.80 

77.12 

76.49 

77.36 


84.62 

84.90 



84.54 

84.82 

84.93 

85.29 


89.51 

89.81 



89.21, 

89.74 

90.72 

90.35 

1 95.00 

94.46 

94.68 



94.41 

94.56 

95.26 

94.87 
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Table A-4.—Body Pressure Orifice Locations 
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For the first 149 runs of the first test in the BTWT, pressure readings at these orifices did not always stabilize. 


Pressure tubing used in this model was 1.016 mm (0.040 in.) o.d. Monel with a 
0.1524 mm (0.006 in.) wall thickness. The major channels for wing pressure tubing were 
machined into the surface. The detailed grooves required to route tubing from the 
orifices to these channels were cut by hand. The pressure orifices were installed normal 
to and flush with the local surface. After installation of the pressure tubing, the grooves 
were filled with solder and brought back to contour by hand-filing to match templates 
prepared by numerically controlled machining. 

Quick disconnects were used at the wing-body junction to reduce the time required for 
installing a different wing. Unfortunately, by the time the cambered-twisted wing was 
installed in the test section, one quick-disconnect block had become worn out due to the 
two previous tests and model checkout. The connection did not seal properly and 
measurements at a series of orifices (x/c from 0.125 through 0.600) on the lower surface 
at 0.80 b/2 were not sufficiently accurate to be used. Data values to be used in the 
integration were obtained by linear spanwise interpolation between adjacent sections. 

The tubing for body pressure orifices was run through the hollow center of the model 
body rather than running it in grooves in the outside contour. Tubing from all the 
orifices was routed through the hollow body to the scanivalves located in the body nose. 
Wiring from the scanivalves was routed through the body to the sting. 


The nose portion of the body was removable to provide access to the fifteen 24-position 
scanivalves. Figure A-1 shows the aft body location of the strain gages that were used to 
measure normal force and pitching moment. 

PRESSURE INSTRUMENTATION 

The model was instrumented with fifteen 24-position scanivalves. Each scani valve 
contained a 103.42-kN/m2 (15-psi) differential Statham, variable resistance, unbonded 
strain gage transducer. These transducers are calibrated against a high accuracy 
standard and, if placed in a temperature-controlled environment, will read within an 
accuracy of 0.1 percent of full scale. The transducers were located inside the model and 
subjected to large temperature excursions. During testing in the Boeing Transonic Wind 
Tunnel (BTWT), temperatures recorded at the scanivalves indicated that the accuracy of 
the readout was 0.75 percent of full-scale capability based on the calibration data. For 
tests in the 9- by 7-ft supersonic leg of the NASA Ames Unitary Wind Tunnel, the 
accuracy of pressure measurements was better than ±0.3 percent, based on the 
maximum temperature measured in the test section. 

During the first test in the BTWT (NASl-12875), two problems were encountered. For 
the first 149 runs, the data filter for one of the scanivalves was inadvertently set at too 
low of a cutoff frequency. This caused a lag that affected five body pressure 
measurements, which produced a maximum error of approximately 0.684 kN/m^ (0.1 psi) 
at an angle of attack of 16° and M = 0.95. Table A-4 identifies the specific data affected. 
During the first half (approximately) of the test, the scanivalve that recorded lower 
surface pressures between the hingelines for the sections at 2y/b = 0.09, 0.20, 0.35, and 
0.50 was intermittent at an angle of attack of 16°. This problem was eventually traced 
to an electrical problem in the strut. Rather than sacrifice all of these data, the 
incorrect measurements were replaced by extrapolating the data from angles of attack 
of 12° and 14°. 
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In the second test in the BTWT (NASl-14962), damage to one of the quick disconnects 
for the wing caused the loss of measurements at a series of orifices (x/c from 0.125 
through 0.600) on the lower surface at 0.80 b/2. Replacement values were obtained by a 
linear-span wise interpolation between the adjacent sections. 

WIND TUNNEL FACILITIES AND CAPABILITIES 

The majority of testing (NASA contracts NASI- 1287 5 and NASl-14962) of this model 
was conducted in the Boeing Transonic Wind Tunnel (BTWT). There was also limited 
testing (NASA contract NASl-14141) of the flat and twisted wings in the 9- by 7-ft 
supersonic leg of the NASA Ames Unitary Wind Tunnel. These two facilities are 
described in some detail. 

BOEING TRANSONIC WIND TUNNEL (BTWT) 

The BTWT is a continuous-flow, closed-circuit, single-return facility with an operating 
range of Mach number from 0.0 to nearly 1.1. The test section is 2.438 by 3.658 by 
4.420 m (8 by 12 by 14,5 ft) with 11.0 percent of the wall area in slots. 

The tunnel layout is shown in figure A-7. The tunnel stagnation pressure is 
atmospheric with a total temperature range of 300 K to 356 K (540® to 640® R). The 
variation with Mach number of Reynolds number based on the mean aerodynamic chord 
(M.A.C.) of this model is shown in figure A-8., which also shows the variation of 
dynamic pressure with Mach number. The 26 856-kW (36 000-hp) wound-rotor induction 
motor in tandem with a 13 428-kW (18 000-hp) synchronous motor provides the power to 
drive a 7.315-m (24-ft) diameter fan up to a maximum speed of 470 rpm. The fan is 
made up of a 5.486-m (18-ft) diameter hub with 72 fixed-pitch fiberglass blades 0.914 m 
(36 in.) long in two stages and directs circuit air through two stages of 67 hollow steel 
stators. 

Data System 

The BTWT data system provides the capabilities of real-time test data acquisition, 
feedback control computation, and display. The data system consists of an Astrodata 
acquisition subsystem and a computing subsystem that uses a Xerox data system 
(XDS 9300) digital computer. The Astrodata system acquires signals from the sensors, 
conditions them, and passes them directly to the computer. Test data (averaged from as 
many as 256 samples per test point) are recorded on a rapid-access data drum. As final 
computations are performed, selected on-line displays are provided on analog X-Y 
plotters and teletypewriters. Real-time computations and displays are performed every 
200 milliseconds for control and test monitoring functions. Any test data may be 
retrieved from rapid-access drum storage and displayed on an oscilloscope. On-line 
programs also provide for the preparation of magnetic tapes for plotting or interfacing 
with off-line programs. Figure A-9 is a schematic of the data acquisition and reduction 
system. 

Mach Number 

Mach number in the BTWT is referenced to the horizontal and lateral center of the test 
section at tunnel station 1000, which was the pitch point of this model (40-percent 
M.A.C.). 
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The pressures used in determining the Mach number, Pg, and p^ are measured through 
permanently positioned sensors. Static pressure Pg is measured by a 103.42-kN/m2 
(15-psi) absolute transducer. A 103.42-kN/m2 (15-psi) differential transducer is used to 
obtain total pressure by measuring (PfPg). These transducers are temperature 
compensated in addition to being in a ±1.11® C (±2®F) environment. Transducer 
performance is checked periodically, and both the static and differential transducers 
have shown a maximum deviation of ±0.02 percent of full scale. 

The static pressure tap is located out of the test section above the ceiling in the 
pressure cap plenum. A correction is made to adjust this static pressure reading to the 
measured test-section-centerline static pressure determined during calibrations at 
station 1000. The tunnel total pressure is obtained from a total pressure probe mounted 
near the tunnel ceiling in the bellmouth throat (fig. A-7). 

Signals from the pressure sensors are fed to the XDS 9300 computer. The XDS system 
computes and updates the Mach display five times per second. Accounting for the entire 
system, calculated Mach number is accurate within ±0.002. Data are recorded only 
when the tunnel is within a preselected Mach tolerance. For this test, a tolerance of 
±0.003 was used. 

Dynamic Pressure 

The dynamic pressure q is computed from the Mach number and the corrected static 
pressure. The estimated tolerance on dynamic pressure is ±95.8 N/m^ (±2.0 psf). 

Angle of Attack 

The angle of attack of the reference point (0.25 M.A.C. for this model) for a 
sting-mounted model is a combination of the input angle measured at the base of the 
sting and several incremental corrections. The input angle of attack is determined by an 
encoder mounted in the strut. This angle is accurate within ±0.02®. This angle is then 
modified by the effects of sting deflection, up-flow, and wall corrections. 

Sting deflections due to load were determined during the calibrations of the strain 
gages, which are mounted on the integral sting body of the model. These deflections are 
known within ±0.02®. The corrections for sting deflection are based on the normal force 
and pitching moment loads obtained during wind-on data acquisition. The sting 
deflection was taken into account when setting test angles of attack, to minimize the 
variation in final angle of attack for the various model configurations. The strain gages 
attached to the sting body of this model have an estimated accuracy of ±5 percent of 
full-scale reading. This means that the sting deflections based on maximum model loads 
were known within ±0.11®. 

During run 55 of the second test in the BTWT, the wiring for the pitching moment gage 
broke, affecting both the normal force and pitching moment measurements and, 
therefore, the calculation of sting deflection under load. For the remainder of that test, 
the model angle was set using the angle of attack as determined by the encoder for the 
most similar previously run configuration. After the test, the normal force and pitching 
moment obtained by integrating the pressure data were used to correct the final angle 
of attack. To verify this procedure, a comparison of these two methods was made using 
data obtained prior to run 55; the results matched within ±0.01®. 
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Up-flow corrections were made based on data obtained from upright and inverted runs 
on a calibration model of similar span. These corrections were less than 0.2®. It is 
generally accepted that the up-flow values are known within ±0.05®. 

A correction to model angle was made for the effect of lift interference for 11-percent 
slotted walls. The lift interference is a function of the ratio of model-to-test section size, 
test section shape, Cj^, and wall geometry. For Cjj = 1.0, this correction is on the order 
of —0.48®. Due to the limited amount of experimental substantiation, the wall correction 
could be in error by ±20 percent. 

NASA AMES UNITARY WIND TUNNEL, 9- BY 7-FT SUPERSONIC LEG 

The 9- by 7-ft supersonic leg of the NASA Ames Unitary Wind Tunnel is a 
continuous-flow, closed-return, variable-density facility with an operating range of 
Mach number from 1.54 to 2.50. (A schematic is shown in fig. A-10.) The tunnel is 
equipped with an asymmetrical sliding-block nozzle and a flexible upper plate; variation 
of the test section Mach number is achieved by translating, in the streamwise direction, 
the fixed-contour block that forms the floor of the nozzle. For this test, the Reynolds 
number was selected as 8.65 by 10® based on the mean aerodynamic chord of this 
model. The test section is 2.74 by 2.13 by 5.49 m (9 by 7 by 18 ft). 

The tunnel air is driven by an 11-stage axial-flow compressor that is powered by four 
variable-speed, wound-rotor induction motors with a combined output of 134 280 kW 
(180 000 hp). Four 850 m® (30 000 ft®) spherical storage tanks provide dry air for tunnel 
pressurization. The temperature is controlled by aftercooling. 

Data System 

The data acquisition system is comprised of a Beckman 210 analog-digital recorder and 
a minicomputer. Output from the Beckman 210 is converted to an acceptable format and 
transmitted by the minicomputer to an IBM 360 computer, which is located in the 
AMES Research Center central computer facility for the processing and preparation of 
final data. This flow is illustrated in figure A-11. 

Angle of Attack 

The angle of attack of the reference point (0.25 M.A.C. for this model) for a 
sting-mounted model is a combination of the input angle at the base of the sting and an 
increment due to sting deflection. The input angle of attack at the base of the sting is 
accurate within 0.02®. 

Sting deflections due to load were determined during the calibration of the strain gages 
mounted on the integral sting body of the model. The corrections for sting deflection are 
based on the normal force and pitching moment loads obtained during wind-on data 
acquisition. The sting deflection was taken into account when setting test angles of 
attack to minimize the variation in final angle of attack for the various model 
configurations. Only a crude calibration of the normal force and pitching moment gages 
was obtained since the force and pitching moment measurements were used primarily 
for calculating sting deflection. Comparison with the integrated pressure results 
indicates that both force and moment measurements may be about 10 percent low, 
which could yield a maximum error in final angle of attack of 0.1®. 
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TESTS AND DATA ACQUISITION 


BOEING TRANSONIC WIND TUNNEL (BTWT) 

Tests 

As previously stated, tests were conducted in the BTWT under two NASA contracts. 
Table A-5 shows the 54 configurations that were tested under contract NAS 1-12875. 
The 12 configurations tested under contract NASI- 14962 are shown in table A-6. Two of 
the configurations were included in both tests to ensure that data from the two tests are 
consistent. Photographs of some of the configurations are shown in figures A- 12 
through A-15; a diagram of the model installation in the BTWT is shown in figure A-16. 

Pressure and total force data were obtained at Mach numbers of 0.40, 0.85, 0.95, and 
1.05 for all configurations and at Mach numbers of 0.70, 1.00, and 1.11 for selected 
configurations. Test angles of attack were from —8® to +16® in 2® increments, however, 
not all angles of attack were included for all configurations and/or Mach numbers. 
Tables A-5 and A-6 show the run numbers for each Mach number and configuration for 
which these data were obtained. 

During the first test, wingtip deflection pictures were taken for representative 
configurations at three Mach numbers to evaluate the stiffness of the wing. These were 
compared to wind-off reference pictures to determine the relative deflection and twist. 
Configurations included the flat and twisted wings, and trailing-edge control surfaces 
deflected +30.2®, 0.0®, and —17.7®. Whereas the tip did deflect (less than 2 cm), the 
change in incidence was negligible even at M = 1.05; and, due to model flexibility, no 
corrections to the data were required. 

Data Repeatability 

Comparisons of data from both tests (NASA CR-165701), show that the data are within 
the tolerances expected for repeat runs during a single test. Therefore, data from both 
tests are compared without regard to the test in which the data were obtained. 

Data Acquisition and Initial Processing 

The pressure data were recorded through the use of fifteen 24-position scanivalves 
located in the fore body of the model. Pressure transducers in the scanivalves measured 
the differential pressure between the local surface pressures and tunnel total pressure. 
Signals from the scanivalves, force and moment data, tunnel parameters, and model 
attitude angle were recorded on the Astrodata system and reduced using the XDS 9300 
computer. 

NASA AMES UNITARY WIND TUNNEL, 9- BY 7-FT SUPERSONIC LEG 
Tests 

Table A- 7 lists the 13 configurations that were tested. Photographs of two of these are 
shown in figures A-17 and A-18; a diagram of the model installation in the test section 
is shown in figure A-19. Pressure and total force data were obtained at Mach numbers 
of 1.70, 2.10 and 2.50 for all configurations. Table A-7 shows the run numbers for each 
Mach number and configuration for which these data were obtained. Test angles of 
attack were from -8® to +14® in 2® increments and +15®. 
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(NASA Contract NASI -12875) 









































































































































Table A-6.— Summary of Subsonic/Transonic Test Conditions by Run Number 
(NASA Contract NASI -14962} 


Mach 

number 

Trailing-edge deflection, degrees 

Full span 

Outboard (inbd=0.0) 

Inboard (outbd=0.0) 

0.0 

8.3 

8.3 

8.3 

Twisted wing 

0.40 

15 

25 

30 

20 

0.70 

14 

24 

29 

19 

0.85 

13 

23 

28 

18 

0.95 

12 

22 

27 

17 

1.00 

11 




1.05 

10 

21 

26 



16 

Cambered-twisted wing, fin off 

0.40 

43 

65 

80 

57 

0.70 

41 




0.85 

40 

62 

78 

59 

0.95 

39 

64 

79 

58 

1.00 

38 




1.05 

37 



61 

77 

55 


Cambered-twisted wing, fin on 


0.40 

49 

70 

75 

54 

0.70 

45 




0.85 

48 

68 

73 

52 

0.95 

47 

69 

74 

53 

1.00 

46 




1.05 

44 

67 

72 

51 
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Data Acquisition and Initial Processing 

The pressure data were recorded through the use of fifteen 24-position scanivalves 
located in the fore body of the model. Pressure transducers in the scanivalves measured 
the differential pressure between the local surface pressures and a known reference 
pressure. Signals from the scanivalves, force and moment data, tunnel parameters, and 
model attitude angle were recorded on the Beckman 210 analog-digital recorder and 
reduced by the Ames staff. 

TRIP STRIP 

A trip strip of no. 60 carborundum grit was used throughout the tests with the 
exception of one series. On the body, the trip strip was 0.32 cm (0.125 in.) wide and 
placed 2.54 cm (1 in.) from the nose. On the wing, it was 0.32 cm (0.125 in.) wide from 
the side of body to the midspan control surface break (0.57 b/2), and tapered to 0.16 cm 
(0.0625 in.) wide at the wingtip. On the upper surface of the wing, the trip strip was 
placed at 15-percent chord; and, on the lower surface, it was placed just aft of the 
location of the leading-edge control surface brackets on the -flat wing (see fig. A-6). 
Density of the grit was 4 to 5 grains per quarter-inch (6 to 8 grains per cm) of trip strip 
length. 

FINAL DATA 

Final data (pressure coefficients, tunnel parameters, and model attitude) were merged 
on magnetic tapes with appropriate configuration and test point identification for 
integration and plotting of these data. 

A detailed description of the data editing and integration procedure are included in 
appendix B. 
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Figure A-1. -General Arrangement and Characteristics 












Figure A-2.—Spanwfse Twist Distribution for the Model Wing 


1.0 



(a) Definition of k. Factor on Basic Camber 


Tip, 2y/b = 1 .00 



Root, 2y/b = 0.0 


Wing reference 








5.994 ( 2 . 360 ) 




All dimensions in centimeters (inches) 



Figure A-4.—Fin Geometry, Section at 0.725 Semispan 
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Figure A-5.-Pressure Orifice Locations 





(a) Schematic 


Tunnel station 



(b) Test Section 

Figure A-7.— Boeing Transonic Wind Tunnel 
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Figure A-9.— Data Acquisition and Reduction System— Boeing Transonic Wind Tunnel 
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mes Unitary Wind Tunnel 






Figure A- 1 1,—Data Acquisition and Reduction System— 9- by 7-ft Supersonic Leg 
of NASA Ames Unitary Wind Tunnel 
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Figure A- 12. -Model in Boeing Transonic Wind Tunnel- Fiat Wing; LE. Deflection, Full 

Span = 0.0°; T.E. Deflection, Full Span = 0.0° (NASA Contract NASI -12875) 
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Figure A-13.—Model in Boeing Transonic Wind Tunnei—Twisted Wing; T.E. Deflection, Full 
Span = 0.0°; (NASA Contract NASI- 14962) 


51 





P inn re A- 14 —Mode! in Boeing Transonic Wind Tunnel— Cam bered-Twisted Wing, Fin Off; T.E 
Deflection, Full Span = O.Cf; (NASA Contract NASI- 14962) 
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Fiaure A-15.-~Model in Boeing Transonic Wind Tunnel— Cambered-Twisted Wing, Fin On; 
T. E. Deflection, Full Span - 0,0° (NASA Contract NASI- 14962) 
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Tunnel ceiling 
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Figure A-16.-Model Installation in Boeing Transonic Wind Tunnel (NASA Contracts NASI -12875 and NASI -14962) 



Figure A-17.—Mode! in 9- by 7-ft Supersonic Leg of NASA Ames Unitary Wind Tunnel— Fiat 
Wing, Rounded L. E. (NASA Contract NASI- 74 14 1) 



Figure A- 18.— Model in 9- by 7-ft Supersonic Leg of NASA Ames Unitary Wind Tunnel- 
Twisted Wing (NASA Contract NASI- 14 14 1) 
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Tunnel (NASA Contract NASI -14141) 



APPENDIX B 


DATA REDUCTION 


DATA EDITING 


There were some cases encountered with these data where the methods of data editing 
available within the integration programs were not adequate. Because the plotting 
program assumes that the geometry is the same for all configurations, and the 
chordwise location of orifices on the various model parts was not absolutely identical, 
points were added as required. Therefore, some interpolations or extrapolations using 
selected orifices were done before the integration program was used. The row of orifices 
on the body at the wing-body intersection was extended in front of the wing and aft of 
the wing by interpolating between the orifices located at 90® and 135®. 


Some specific problems with the data acquisition systems required the replacement of 
some data for parts of tests. These are identified in appendix A, 

Several methods were introduced into the integration program to replace or add data 
points to account for: 

- Plugged or leaking orifices, or bad data points 
Extrapolating the data to leading and trailing edges 

- Hingeline discontinuities in the pressure data 


These procedures were selected by code for each point. The codes are described in the 
following list and are illustrated in figure B-1. An additional use of these codes is to 
ensure that only measured pressure data (CODE^ = 0) are identified with symbols on 
the plots. The subscript i identifies the position of the point from the leading edge of the 
upper or lower surface of the section. 

IF CODEj = 0, use pressure as entered on tape (measured pressure) 

= 20, use as entered on tape (previously replaced value) 


= 1, interpolate from adjacent points 
= 2, extrapolate from two preceding points 
= 3, extrapolate from two following points 
= 4, set equal to preceding point 
= 5, set equal to following point 
= 6, interpolate using points (i-2) and (i + 1) 
= 7, interpolate using points (i-1) and (i+2) 
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IF CODEj = negative of above, evaluate as above but average with corresponding point 
on opposite surface; used for leading and trailing edges of section only 

Editing of the pressure data is done in the following order: 

1. Each section is done separately. 

2. Each surface (upper or lower) per section is done in the following sequence; 

a. Starting at leading edge, points with codes of 1, 2, and 4. 

b. Starting at trailing edge, points with codes of 3, 5, 6, and 7. 

3. Leading- and trai ling-edge points with negative codes are evaluated. Both upper 
and lower surface codes need not be negative and need not be the same negative 
code. 


CALCULATION OF NET PRESSURE COEFFICIENTS 
The net lift distribution on the section is calculated by: 


C = C C 

p,net p, lower ~ ^p, upper 


(B-1) 


INTEGRATION OF PRESSURE DATA 


To account for the effects on integrated coefficients of the deflected control surfaces, 
each streamwise section (of which there are NSECT) is divided into segpnents (of which 
there are NSEG). These segments are the leading-edge control surface, wing box, and 
trailing-edge control surface. The upper and lower surfaces of each are integrated 
separately over the number of points available ((number of orifices + 2) = NPl), and 
are based on the segment chord length c. Sign conventions are shown in the following 
sketch. The equations, which use a rectangular integration process, follow. 
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SEGMENT COEFFICIENTS 

Integration of the pressures for each segment per surface per section is the first step. 


o Normal force coefficient 


NPl 


S hvwMfe-ei,] 

^n,s,net ^ ^n,s, lower “ ^n,s, upper 


o Chord force coefficient C 


c,s 


,s = 0.5 E[(C,,)^.(Cp)^_||>-(l),.,] 


Cc a — 


^c,s,net ^c,s, upper “ ^c,s, lower 

o Pitching moment coefficient about segment leading edge g 


NP 


Cn,,s = 0.5i: [(Cp),.(Cp)._,; 


( tUtX 




L ^ i-1 


2.0 


[©r©,J 


NP 


- o,«E [(s ),.( c,)^j 


^m,s,net ~ ^m,s, upper g lower 


0 Pitching moment coefficient about 0.25 c of segment Cm.25c,s 


^m.25c,s “ ^m,s, ^n,s 


(B-2) 

(B-3) 

(B-4) 

(B-5) 

(B-6) 

(B-7) 

(B-8) 
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SECTION COEFFICIENTS 


Total section coefficients are obtained by summing the segment coefficients, taking into 
account both the segment deflections as defined in the following sketch and the segment 
chord lengths. These coefficients are based on the section chord length c. 


-8 



This sign convention for leading-edge deflection is used only in summation 
of coefficients. 



o Normal force coefficient C, 


NSEG 


NSEG 


C„=E Gn,s),(^J cosV E 

o Pitching moment coefficient about section leading edge Cjj^ 


(B-9) 


C 


m 


NSEG 


g (t); + (Cn,s) , (‘ 0 - cos S ,) + (c,,,) ^ sin 6, 


NSEG / V 


(B-10) 


where 


Xl.E.,s 

^L.E. 


is segment chord length, cm 
is section chord length, cm 

is deflection of segment relative to section chord plane, leading edge 
up, degrees 

is leading edge of segment, cm 
is leading edge of section, cm 


0 Pitching moment coefficient about 0.25 c of section 25c 

m.25c “ 


(B-11) 
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TOTAL SURFACE COEFFICIENTS 


To obtain total surface coefficients, the assumption is made that the section coefficients 
apply for a finite distance on both sides of each row of orifices. The equations for total 
surface coefficients are as follows: 


o Normal force coefficient Cj^ 


NSECT 

Cn=J 


o Bending moment coefficient Cg 

NSECT 

^B=-S(b7r) 


(B-12) 


(B-13) 


o Pitching moment coefficient about 0.25 M.A.C. 


Cm-§F 


NSECT , 

,5, ('^m)k(V)k Cn)k(Sh)k [’‘ref-(>^L.E.),,] 


(B-14) 


where 


c is reference chord for pitching moment, cm 

Xj.gf is reference station for pitching moment, cm (0.25 M.A.C.) 

^L.E. i® leading edge of section chord, cm 

b/2 is reference length for bending moment, cm 
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DETERMINATION OF GEOMETRIC CONSTANTS REQUIRED FOR 
INTEGRATION 

To obtain total surface coefficients, the assumption is made that the section coefficients 
apply for a finite distance on both sides of each row of orifices. The input geometry 
required to calculate the areas, and products of area and length required for the 
summation of total surface coefficients, is shown in the following sketch. 



0 Section area: 


X'p + tan Ay £ y 


/* yout /“ 

I I dy dx 

Ain AL + tanAL,E.y 


(*■■*" ^T.E. '^L.E.) (yout^ '^ 111 ^) (B-15) 


o Product of section area and mean chord: 

•Vout /•’‘T + '^T.E y 


S|,c = 


/ j'oiit r 

I c dy dx 

yjn ^L.E.y 


= v(yout-yin) (‘an At.E. - ''l.E.) (yout^*yin^) 
( tan Ay £ - tan Al e.) ' / 3 3 ) 

+ 3 Q vout ^in / 


(B-16) 
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o 


Product of section area and bending moment arm: 


/-Vout + 

s^y = / / (y-Vrefrydx 

^ c,.-(tan A-t e - tan A|^ g )yref . y,^2) 

o Product of section area and leading-edge coordinate: 


/•yout r XT + tanATjgy 

ShX =/ / ’^L.E.dydx 

■'yin ■'^’‘L + ^'L.E.y 


= xlc 


+ tan A 


tan Ai^ E c^. + (tan Ay ^ - tan Al e. )/.. 2 

'r(^out“yin/ 2.0 ^ 


Vi 


in 


( tan A-t £ - tan ^ 3^3^ 

L.E. 3.0 Vyoiit ^in / 


o Total surface reference area: 


NSECT 


S = 


?. (^h)i 


k=] 


M.A.C. and X coordinate of M.A.C. leading edge: 

NSECT 

__ 1 


“""s £ (V) 


k=l 


NSECT 

■L.E.,M.A.C.=-? L (®hx)| 


k=l 


The required integration constants for the wing and body are shown in table B-1. 


(B-17) 

(B-18) 

(B-19) 

(B-20) 

(B-21) 
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Table B-1.— Integration Constants 


Reference area = 3128.45 cm^ 
M.A.C. = 75.311 cm 

Half span = 50.80 cm 


Pitching moment referenced to 0.25 M.A.C. 

Bending moment referenced to 0.086 = 4.374 cm) 


Wing 


LE. of M.A.C. at B.S. 87.760 cm 


2v/b 

Ay 

Area 

Area • chord 

Area • (y-Vref ) 

(b/2) 

cm2 

cm^ 

cm2 

0.09 

0.0425 

219.69 

22 357. 

167. 

0.20 

0.1575 

733.51 

67 415. 

4 206. 

0.35 

0.1500 

580.54 

44 374. 

7 857. 

0.50 

0.1400 

437.93 

27 084. 

9 148. 

0.65 

0.1600 

377.64 

17 722. 

10 729. 

0.80 

0.1300 

210.35 

6 794. 

7 528. 

0.93 

0.1400 

129.79 

2 487. 

5 505. 


Body 


Longitudinal 

Area 

Area • L 

section 

cm2 

cm2 

1 

356.61 

81 258. 

2 

504.32 

114916. 

3 

- - 

70.94 

16 164. 
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Extrapolated 





\ 



code = -2 
(similar 
at leading 
edge when 
code = —3) 


Figure B-1.— Codes Used to Interpolate and Extrapolate 
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Figure 1 .—General Arrangement and Characteristics 
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Figure 3. -Cambered-T wisted Wing Section Geometry 
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Figure 4.— Pressure Orifice Locations 


Figure S.-Surface Paneling for Arrow Wing 



Fraction of local chord - x/c 


Figure 6.— Comparison of Nondimensionai Definition of the Trailing Edge of the Airfoil, Fiat Wing 



Upper surface pressure coefficient Upper surface pressure coefficient 


Fraction of semispan .09 



Fraction of semispan .20 



Fraction of semispan .65 



Fraction of semispan .80 



Figure 7.— Comparison of Predicted Pressure Distributions for Two Airfoil Definitions, Fiat 
Wing, Rounded L.E., M = 0.40, ol= 16° 
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upper surface pressure coefficient Upper surface pressure coefficient 
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Fraction of semispan .50 



Fraction of semispan .93 



Figure 7.— (Concluded) 
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Upper surface pressure coefficient Upper surface pressure coefficient 


Fraction of semispan .09 



Fraction of semispan .65 



Figure 9.— Effect of Initial Core Radius in White's 


Fraction of semispan .20 



Fraction of local chord 


Fraction of semispan .80 



Fraction of local chord 


Fiat Wing, Rounded L.E., M = 0.40, a = 16' 
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Figure 9.—( Concluded) 
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Figure 10.~f Concluded) 
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Upper surface pressure coefficient Upper surface pressure coefficient 


Fraction of semispan .09 



Fraction of semispan .20 



Fraction of semispan .65 



Fraction of semispan .80 



Figure 1 1 .-Theory-to-Experiment Comparison of Wing Upper Surface Pressure 
Distributions— Flat Wing, Rounded L.E., M = 0.40, a = 16F 
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Fraction of semispan .50 








Upper surface pressure coefficient Upper surface pressure coefficient 


Fraction of semispan .09 Fraction of semispan .20 



Fraction of semispan .65 



Fraction of semispan .80 



Figure 1 2.— Theory-to- Experiment Comparison of Wing Upper Surface Pressure 
Distributions— Twisted Wing, Rounded L.E., M = 0.40, a = 16° 
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Fraction of semispan .35 
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Figure 12.~( Concluded) 
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Upper surface pressure coefficient Upper surface pressure coefficient 


Fraction of semispan .09 Fraction of semispan .20 



Fraction of local chord Fraction of local chord 


Fraction of semispan .65 Fraction of semispan .80 



Fraction of local chord Fraction of local chord 


Figure 1 3.— Theory-to- Experiment Comparison of Wing Upper Surface Pressure 

Distributions— Cambered-Twisted Wing, Rounded L.E., M = 0.40, a = 16° 
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Fraction of semispan .35 


Fraction of semispan .50 







Fraction of semispan .65 Fraction of semispan .80 



Fraction of lo.cal chord Fraction of local chord 


Figure 1 4. -Theory-to-Experiment Comparison of Wing Upper Surface Pressure 
Distributions— Fiat Wing, Rounded L.E., M = 0.40, a = 8° 
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were evaluated under this contract. These three tasks are documented separately. 


This report describes the evaluation of R. P. White’s (RASA Division of Systems 
Research Laboratories) separated-flow method (Task I). This method was developed for 
moderately swept wings with multiple, constant-strength vortex systems. The flow on 
the highly swept wing used in this evaluation is characterized by a single vortex system 
of continuously varying strength. The data comparisons show, as currently formulated, 
that this method does not predict the pressure distribution on this highly swept wing. 
The evaluation of The Boeing Company’s Three-Dimensional Leading-Edge Vortex 
(LEV) code (Task III) is presented in NASA CR-3642. The development and evaluation 
of a semi-empirical method to predict pressure distributions on a deformed wing from an 
experimental data base (Task II) is described in NASA CR-3641. 

These evaluations were based on a large experimental data base (for three wing 
shapes), which was obtained under previous NASA contracts: NASl-12875, 
NASl-14141, and NASl-14962. Linear theoretical methods were also evaluated under 
these contracts. 
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